Eradicating typological thinking in prokaryotic systematics and evolution.

In 1982, John Maynard Smith called for an evolutionary "New Synthesis" specific for prokaryotes, observing that "population thinking has been well developed for fully half a century, but has yet to be adopted by microbiology" (Maynard Smith 1982). Twenty-seven years later, typological thinking (population thinking's antithesis) still dominates the field. Evidence for this includes the continuing debates on the reality of prokaryotic species, the value of the term "prokaryote," and the significance of the tree of life (TOL). In each case, the unexpected prevalence of interlineage transfer of genetic information has been (or should now be) the catalyst for the final Darwinization of our discipline. With examples from phylogenomics, I argue that "species," "domains," and the "TOL" are reifications that we can do without, especially as genomics dissolves into metagenomics.

Integron-associated gene cassettes in Halifax Harbour: assessment of a mobile gene pool in marine sediments.

The integron/gene cassette systems identified in bacteria comprise a class of genetic elements that allow adaptation by acquisition of gene cassettes. Integron gene cassettes have been shown to facilitate the spread of drug resistance in human pathogens but their role outside a clinical setting has not been explored extensively. We sequenced 2145 integron gene cassettes from four marine sediment samples taken from the vicinity of Halifax Nova Scotia, Canada, increasing the number of gene cassettes obtained from environmental microbial communities by 10-fold. Sequence analyses reveals that the majority of these cassettes encode novel proteins and that this study is consistent with previous claims of high cassette diversity as we estimate a Chao1 diversity index of approximately 3000 cassettes from these samples. The functional distribution of environmental cassettes recovered in this study, when compared with that of cassettes from the only other source with significant sampling (Vibrio genomes) suggests that alternate selection regimes might be acting on these two gene pools. The majority of cassettes recovered in this study encode novel, unknown proteins. In instances where we obtained multiple alleles of a novel protein we demonstrate that non-synonymous versus synonymous substitution rates ratios suggest relaxed selection. Cassette-encoded proteins with known homologues represent a variety of functions and prevalent among these are isochorismatases; proteins involved in iron scavenging. Phylogenetic analysis of these isochorismatases as well as of cassette-encoded acetyltransferases reveals a patchy distribution, suggesting multiple sources for the origin of these cassettes. Finally, the two most environmentally similar sample sites considered in this study display the greatest overlap of cassette types, consistent with the hypothesis that cassette genes encode adaptive proteins.

Alternative methods for concatenation of core genes indicate a lack of resolution in deep nodes of the prokaryotic phylogeny.

It has recently been proposed that a well-resolved Tree of Life can be achieved through concatenation of shared genes. There are, however, several difficulties with such an approach, especially in the prokaryotic part of this tree. We tackled some of them using a new combination of maximum likelihood-based methods, developed in order to practice as safe and careful concatenations as possible. First, we used the application concaterpillar on carefully aligned core genes. This application uses a hierarchical likelihood-ratio test framework to assess both the topological congruence between gene phylogenies (i.e., whether different genes share the same evolutionary history) and branch-length congruence (i.e., whether genes that share the same history share the same pattern of relative evolutionary rates). We thus tested if these core genes can be concatenated or should be instead categorized into different incongruent sets. Second, we developed a heat map approach studying the evolution of the phylogenetic support for different bipartitions, when the number of sites of different phylogenetic quality in the concatenation increases. These heatmaps allow us to follow which phylogenetic signals increase or decrease as the concatenation progresses and to detect emerging artifactual groupings, that is, groups that are more and more supported when more and more homoplasic sites are thrown in the analysis. We showed that, as far as 7 major prokaryotic lineages are concerned, only 22 core genes can be said to be congruent and can be safely concatenated. This number is even smaller than the number of genes retained to reconstruct a "Tree of One Per Cent." Furthermore, the concatenation of these 22 markers leads to an unresolved tree as the only groupings in the concatenation tree seem to reflect emerging artifacts. Using concatenated core genes as a valid framework to classify uncharacterized environmental sequences can thus be misleading.

Visualizing and assessing phylogenetic congruence of core gene sets: a case study of the gamma-proteobacteria.

Here, we address a much-debated topic: is there or is there not an organismal tree of gamma-proteobacteria that can be unambiguously inferred from a core of shared genes? We apply several recently developed analytical methods to this problem, for the first time. Our heat map analyses of P values and of bootstrap bipartitions show the presence of conflicting phylogenetic signals among these core genes. Our synthesis reconstruction suggests that at least 10% of these genes have been laterally transferred during the divergence of the gamma-proteobacteria, and that for most of the rest, there is too little phylogenetic signal to permit firm conclusions about the mode of inheritance. Although there is clearly a central tendency in this data set (it is far from random), lateral gene transfers cannot be ruled out. Instead of an organismal tree, we propose that these core genes could be used to define a more subtle and partially reticulated pattern of relationships.

The genome of Salinibacter ruber: convergence and gene exchange among hyperhalophilic bacteria and archaea.

Saturated thalassic brines are among the most physically demanding habitats on Earth: few microbes survive in them. Salinibacter ruber is among these organisms and has been found repeatedly in significant numbers in climax saltern crystallizer communities. The phenotype of this bacterium is remarkably similar to that of the hyperhalophilic Archaea (Haloarchaea). The genome sequence suggests that this resemblance has arisen through convergence at the physiological level (different genes producing similar overall phenotype) and the molecular level (independent mutations yielding similar sequences or structures). Several genes and gene clusters also derive by lateral transfer from (or may have been laterally transferred to) haloarchaea. S. ruber encodes four rhodopsins. One resembles bacterial proteorhodopsins and three are of the haloarchaeal type, previously uncharacterized in a bacterial genome. The impact of these modular adaptive elements on the cell biology and ecology of S. ruber is substantial, affecting salt adaptation, bioenergetics, and photobiology.

Do orthologous gene phylogenies really support tree-thinking?

BACKGROUND: Since Darwin's Origin of Species, reconstructing the Tree of Life has been a goal of evolutionists, and tree-thinking has become a major concept of evolutionary biology. Practically, building the Tree of Life has proven to be tedious. Too few morphological characters are useful for conducting conclusive phylogenetic analyses at the highest taxonomic level. Consequently, molecular sequences (genes, proteins, and genomes) likely constitute the only useful characters for constructing a phylogeny of all life. For this reason, tree-makers expect a lot from gene comparisons. The simultaneous study of the largest number of molecular markers possible is sometimes considered to be one of the best solutions in reconstructing the genealogy of organisms. This conclusion is a direct consequence of tree-thinking: if gene inheritance conforms to a tree-like model of evolution, sampling more of these molecules will provide enough phylogenetic signal to build the Tree of Life. The selection of congruent markers is thus a fundamental step in simultaneous analysis of many genes. RESULTS: Heat map analyses were used to investigate the congruence of orthologues in four datasets (archaeal, bacterial, eukaryotic and alpha-proteobacterial). We conclude that we simply cannot determine if a large portion of the genes have a common history. In addition, none of these datasets can be considered free of lateral gene transfer. CONCLUSION: Our phylogenetic analyses do not support tree-thinking. These results have important conceptual and practical implications. We argue that representations other than a tree should be investigated in this case because a non-critical concatenation of markers could be highly misleading.

How big is the iceberg of which organellar genes in nuclear genomes are but the tip?

As more and more complete bacterial and archaeal genome sequences become available, the role of lateral gene transfer (LGT) in shaping them becomes more and more clear. Over the long term, it may be the dominant force, affecting most genes in most prokaryotes. We review the history of LGT, suggesting reasons why its prevalence and impact were so long dismissed. We discuss various methods purporting to measure the extent of LGT, and evidence for and against the notion that there is a core of never-exchanged genes shared by all genomes, from which we can deduce the "true" organismal tree. We also consider evidence for, and implications of, LGT between prokaryotes and phagocytic eukaryotes.

Reconstructing/deconstructing the earliest eukaryotes: how comparative genomics can help.

We could reconstruct the evolution of eukaryote-specific molecular and cellular machinery if some living eukaryotes retained primitive cellular structures and we knew which eukaryotes these were. It's not clear that either is the case, but the expanding protist genomic database could help us in several ways.

Defining the core of nontransferable prokaryotic genes: the euryarchaeal core.

If lateral gene transfer (LGT) has affected all genes over the course of prokaryotic evolution, reconstruction of organismal phylogeny is compromised. However, if a core of genes is immune to transfer, then the evolutionary history of that core might be our most reliable guide to the evolution of organisms. Such a core should be preferentially included in the subset of genes shared by all organisms, but where universally conserved genes have been analyzed, there is too little phylogenetic signal to allow determination of whether or not they indeed have the same history (Hansmann and Martin 2000; Teichmann and Mitchison 1999). Here we look at a more restricted set, 521 homologous genes (COGs) simultaneously present in four sequenced euryarchaeal genomes. Although there is overall little robust phylogenetic signal in this data set, there is, among well-supported trees, strong representation of all three possible four-taxon topologies. "Informational" genes seem no less subject to LGT than are "operational genes," within the euryarchaeotes. We conclude that (i) even in this collection of conserved genes there has been extensive LGT (orthologous gene replacement) and (ii) the notion that there is a core of nontransferable genes (the "core hypothesis") has not been proven and may be unprovable.

Gene duplication and the evolution of group II chaperonins: implications for structure and function.

Chaperonins are multisubunit protein-folding assemblies. They are composed of two distinct structural classes, which also have a characteristic phylogenetic distribution. Group I chaperonins (called GroEL/cpn60/hsp60) are present in Bacteria and eukaryotic organelles while group II chaperonins are found in Archaea (called the thermosome or TF55) and the cytoplasm of eukaryotes (called CCT or TriC). Gene duplication has been an important force in the evolution of group II chaperonins: Archaea possess one, two, or three homologous chaperonin subunit-encoding genes, and eight distinct CCT gene families (paralogs) have been described in eukaryotes. Phylogenetic analyses indicate that while the duplications in archaeal chaperonin genes have occurred numerous times independently in a lineage-specific fashion, the eight different CCT subunits found in eukaryotes are the products of duplications that occurred early and very likely only once in the evolution of the eukaryotic nuclear genome. Analyses of CCT sequences from diverse eukaryotic species reveal that each of the CCT subunits possesses a suite of invariant subunit-specific amino acid residues ("signatures"). When mapped onto the crystal structure of the archaeal chaperonin from Thermoplasma acidophilum, these signatures are located in the apical, intermediate, and equatorial domains. Regions that were found to be variable in length and/or amino acid sequence were localized primarily to the exterior of the molecule and, significantly, to the extreme tip of the apical domain (the "helical protrusion"). In light of recent biochemical and electron microscopic data describing specific CCT-substrate interactions, our results have implications for the evolution of subunit-specific functions in CCT.

Bacterial origin for the isoprenoid biosynthesis enzyme HMG-CoA reductase of the archaeal orders Thermoplasmatales and Archaeoglobales.

The enzyme 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMG-CoA reductase or HMGR) fulfills an essential role in archaea, as it is required for the synthesis of isoprenoid ethers, the main component of archaeal cell membranes. There are two clearly homologous but structurally different classes of the enzyme, one found mainly in eukaryotes and archaea (class 1), and the other found in bacteria (class 2). This feature facilitated the identification of several cases of interdomain lateral gene transfer (LGT), in particular, the bacterial origin for the HMGR gene from the archaeon Archaeoglobus fulgidus. In order to investigate if this LGT event was recent and limited in its scope or had a broad and long-term impact on the recipient and its related lineages, the HMGR gene was amplified and sequenced from a variety of archaea. The survey covered close relatives of A. fulgidus, the only archaeon known prior to this study to possess a bacterial-like HMGR; representatives of each main euryarchaeal group were also inspected. All culturable members of the archaeal group Archaeoglobales were found to display an HMGR very similar to the enzyme of the bacterium Pseudomonas mevalonii. Surprisingly, two species of the genus Thermoplasma also harbor an HMGR of bacterial origin highly similar to the enzymes found in the Archaeoglobales. Phylogenetic analyses of the HMGR gene and comparisons to reference phylogenies from other genes confirm a common bacterial origin for the HMGRs of Thermoplasmatales and Archaeoglobales. The most likely explanation of these results includes an initial bacteria-to-archaea transfer, followed by a another event between archaea. Their presence in two divergent archaeal lineages suggests an important adaptive role for these laterally transferred genes.

The complete genome of the crenarchaeon Sulfolobus solfataricus P2.

The genome of the crenarchaeon Sulfolobus solfataricus P2 contains 2,992,245 bp on a single chromosome and encodes 2,977 proteins and many RNAs. One-third of the encoded proteins have no detectable homologs in other sequenced genomes. Moreover, 40% appear to be archaeal-specific, and only 12% and 2.3% are shared exclusively with bacteria and eukarya, respectively. The genome shows a high level of plasticity with 200 diverse insertion sequence elements, many putative nonautonomous mobile elements, and evidence of integrase-mediated insertion events. There are also long clusters of regularly spaced tandem repeats. Different transfer systems are used for the uptake of inorganic and organic solutes, and a wealth of intracellular and extracellular proteases, sugar, and sulfur metabolizing enzymes are encoded, as well as enzymes of the central metabolic pathways and motility proteins. The major metabolic electron carrier is not NADH as in bacteria and eukarya but probably ferredoxin. The essential components required for DNA replication, DNA repair and recombination, the cell cycle, transcriptional initiation and translation, but not DNA folding, show a strong eukaryal character with many archaeal-specific features. The results illustrate major differences between crenarchaea and euryarchaea, especially for their DNA replication mechanism and cell cycle processes and their translational apparatus.

Microbial genomes: dealing with diversity.

We have now complete genome sequences of several pairs of closely related prokaryotes (conspecific strains or congeneric species). Surprisingly, even strains of the same species can differ by as much as 20% in gene content. Conceptual and methodological approaches for dealing with such diversity are now being developed, and should transform microbial genomics.

Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria.

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.

Class I release factors in ciliates with variant genetic codes.

In eukaryotes with the universal genetic code a single class I release factor (eRF1) most probably recognizes all stop codons (UAA, UAG and UGA) and is essential for termination of nascent peptide synthesis. It is well established that stop codons have been reassigned to amino acid codons at least three times among ciliates. The codon specificities of ciliate eRF1s must have been modified to accommodate the variant codes. In this study we have amplified, cloned and sequenced eRF1 genes of two hypotrichous ciliates, Oxytricha trifallax (UAA and UAG for Gln) and Euplotes aediculatus (UGA for Cys). We also sequenced/identified three protist and two archaeal class I RF genes to enlarge the database of eRF1/aRF1s with the universal code. Extensive comparisons between universal code eRF1s and those of Oxytricha, Euplotes, and Tetrahymena which represent three lineages that acquired variant codes independently, provide important clues to identify stop codon-binding regions in eRF1. Domain 1 in the five ciliate eRF1s, particularly the TASNIKS heptapeptide and its adjacent region, differs significantly from domain 1 in universal code eRF1s. This observation suggests that domain 1 contains the codon recognition site, but that the mechanism of eRF1 codon recognition may be more complex than proposed by Nakamura et al. or Knight and Landweber.

Tubulins in Trichomonas vaginalis: molecular characterization of alpha-tubulin genes, posttranslational modifications, and homology modeling of the tubulin dimer.

We have isolated and analysed an alpha-tubulin-encoding gene (atub1) in an early-diverging eukaryote, Trichomonas vaginalis. The complete atub1 open reading frame included 1.356 bp encoding a polypeptide of 452 amino-acyl residues. A second alpha-tubulin gene (atub2) was amplified by PCR using primers derived from consensus alpha-tubulin amino acid sequences. Both T. vaginalis alpha-tubulin sequences showed high identity to those described in other parabasalids (94.4%-97.3%), and exhibited a high degree of similarity to sequences from Metazoa (such as pig brain) and diplomonads (such as Giardia). Despite large evolutionary distances previously observed between trichomonads and mammals, the three-dimensional model of the T. vaginalis tubulin dimer was very similar to that of pig brain. Possible correlations between alpha-tubulin sequences and posttranslational modifications (PTMs) were examined. Our observations corroborated previous data obtained in T. vaginalis using specific anti-PTMs antibodies. As described in the related species Tritrichomonas mobilensis, microtubules are likely acetylated, non-tyrosinated, glutamylated, and non-glycylated in T. vaginalis. Evolutionary considerations concerning the time of appearance of these tubulin PTMs are also discussed since trichomonads are potentially one of the earliest diverging eukaryotic lineages.

An updated and comprehensive rRNA phylogeny of (crown) eukaryotes based on rate-calibrated evolutionary distances.

Recent experience with molecular phylogeny has shown that all molecular markers have strengths and weaknesses. Nonetheless, despite several notable discrepancies with phylogenies obtained from protein data, the merits of the small subunit ribosomal RNA (SSU rRNA) as a molecular phylogenetic marker remain indisputable. Over the last 10 to 15 years a massive SSU rRNA database has been gathered, including more then 3000 complete sequences from eukaryotes. This creates a huge computational challenge, which is exacerbated by phenomena such as extensive rate variation among sites in the molecule. A few years ago, a fast phylogenetic method was developed that takes into account among-site rate variation in the estimation of evolutionary distances. This "substitution rate calibration" (SRC) method not only corrects for a major source of artifacts in phylogeny reconstruction but, because it is based on a distance approach, allows comprehensive trees including thousands of sequences to be constructed in a reasonable amount of time. In this study, a nucleotide variability map and a phylogenetic tree were constructed, using the SRC method, based on all available (January 2000) complete SSU rRNA sequences (2551) for species belonging to the so-called eukaryotic crown. The resulting phylogeny constitutes the most complete description of overall eukaryote diversity and relationships to date. Furthermore, branch lengths estimated with the SRC method better reflect the huge differences in evolutionary rates among and within eukaryotic lineages. The ribosomal RNA tree is compared with a recent protein phylogeny obtained from concatenated actin, alpha-tubulin, beta-tubulin, and elongation factor 1-alpha amino acid sequences. A consensus phylogeny of the eukaryotic crown based on currently available molecular data is discussed, as well as specific problems encountered in analyzing sequences when large differences in substitution rate are present, either between different sequences (rate variation among lineages) or between different positions within the same sequence (among-site rate variation).

A kingdom-level phylogeny of eukaryotes based on combined protein data.

Current understanding of the higher order systematics of eukaryotes relies largely on analyses of the small ribosomal subunit RNA (SSU rRNA). Independent testing of these results is still limited. We have combined the sequences of four of the most broadly taxonomically sampled proteins available to create a roughly parallel data set to that of SSU rRNA. The resulting phylogenetic tree shows a number of striking differences from SSU rRNA phylogeny, including strong support for most major groups and several major supergroups.

Origin and evolution of eukaryotic chaperonins: phylogenetic evidence for ancient duplications in CCT genes.

Chaperonins are oligomeric protein-folding complexes which are divided into two distantly related structural classes. Group I chaperonins (called GroEL/cpn60/hsp60) are found in bacteria and eukaryotic organelles, while group II chaperonins are present in archaea and the cytoplasm of eukaryotes (called CCT/TriC). While archaea possess one to three chaperonin subunit-encoding genes, eight distinct CCT gene families (paralogs) have been characterized in eukaryotes. We are interested in determining when during eukaryotic evolution the multiple gene duplications producing the CCT subunits occurred. We describe the sequence and phylogenetic analysis of five CCT genes from TRICHOMONAS: vaginalis and seven from GIARDIA: lamblia, representatives of amitochondriate protist lineages thought to have diverged early from other eukaryotes. Our data show that the gene duplications producing the eight CCT paralogs took place prior to the organismal divergence of TRICHOMONAS: and GIARDIA: from other eukaryotes. Thus, these divergent protists likely possess completely hetero-oligomeric CCT complexes like those in yeast and mammalian cells. No close phylogenetic relationship between the archaeal chaperonins and specific CCT subunits was observed, suggesting that none of the CCT gene duplications predate the divergence of archaea and eukaryotes. The duplications producing the CCTdelta and CCTepsilon subunits, as well as CCTalpha, CCTbeta, and CCTeta, are the most recent in the CCT gene family. Our analyses show significant differences in the rates of evolution of archaeal chaperonins compared with the eukaryotic CCTs, as well as among the different CCT subunits themselves. We discuss these results in light of current views on the origin, evolution, and function of CCT complexes.